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Abstract 

We conduct molecular dynamics simulations of 887 and 1389-atom decahedral platinum nanopar- 
ticles using an embedded atom potential. By constructing micro canonical caloric curves, we identify 
structural transitions from decahedral to fee in the particles prior to melting. The transitions take 
place during phase coexistence and appear to occur via melting of the decahedral structure and 
subsequent recrystallisation into the fee structure. 



I. INTRODUCTION 



It has been known for some time that nanoscale particles of fee metals can adopt non- 
crystalline structures, such as icosahedra and truncated decahedra, as they trade internal 
strain for energetically favourable close-packed surface facets^. Generally it is found that 
icosahedra are favoured at small sizes and truncated decahedra at intermediate sizes before 
the particles eventually adopt their crystalline fee structures^. However, much less is known 
about the thermal stability of such non-crystalline phases. Molecular dynamics simulations 
of Lennard- Jones clusters have found that the non-crystalline phases became more stable at 
temperatures approaching the melting point i.e. the size threshold at which these clusters 
prefered the non-crystalline phases tended to increase^. However recent experiments by 
Koga et al 5 where free gold clusters were annealed near the melting point found that small 
particles underwent a transformation from icosahedral to decahedral structures, a finding 
inconsistent with earlier theoretical considerations. 

Other simulations have suggested that solid-liquid phase coexistence may play a role in 
structural transitions 6,7 . There are now several instances where structural transitions have 
been seen in partially melted clusters. In one instance, a transition from an icosahedral 
to decahedral structure, driven by a preference for the melt to wet the (lOO)-facets of the 
truncated decahedron, was observed in a 1415-atom nickel cluster—. In another instance, a 
transition from decahedral structures to icosahedral structures was seen in partially melted 
palladium clusters 7 . The underlying cause of this transition was not determined but it 
was suggested that it may have been due to the preference of recrystallisation kinetics for 
icosahedral structures after fluctuations in the liquid fraction of the particles. These studies 
suggest that phase coexistence (and possibly surface melting) may play an important role 
in structural transitions in nanoparticles. 

HRTEM studies of platinum nanoparticles have found that they showed a preference for 
fee structures down to quite small sizeA Yet molecular dynamics studies using empirical 
potentials^ suggest that non-crystalline structures ought to be stable at low temperatures 
at least up to sizes of several thousand atoms. This inconsistency between theory and 
experiment may be due either to the favouring of fee structures by growth kinetics 10 or to 
the thermal instability of non-crystalline structures^. Platinum clusters are of particular 
interest because of their applications to catalysi o 12113 . Catalytic activity depends critically 
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on the surface structure of the nanoparticle and thus the thermal stability of small fee and 
non- crystalline particles is directly relevant to their catalytic efficiency 

In this paper, we report on molecular dynamics simulations of platinum clusters using 
an embedded atom method (EAM) potential^ where we observe structural transitions from 
non- crystalline truncated decahedral structures to fee stuctures prior to melting. As in the 
case of the transitions seen in Pd clusters^, the structural stransformations appear to occur 
via a melting and recrystallisation process although in this case the kinetics are an order of 
magnitude slower than in Pd. 

II. METHODOLOGY 

We use molecular dynamics to simulate the caloric curves of Pt clusters using an EAM 
potential^. We attempted to identify the stable structures at zero temperature by comparing 
the potential energies of closed-shell truncated octahedra (including variants such as the 
(TO) - and (TO) + structures^), cuboctahedra, Mackay icosahedra, and Marks and Ino 
decahedral. For each structural sequence the potential energy per atom, E, was fitted to 
the expression E(N) = A + BN^ 1 ^ 3 + CN~ 2 ^ 3 . We find that Marks decahedra are favoured 
at relatively small sizes (887 and 1389-atoms), with fee structures being favoured at sizes 
above this. At intermediate sizes the most stable fee structure is the (TO) + although at 
large sizes the truncated octahedra appear to be favoured. In Fig. ^ we only show the 
relaxed energies and fits for the the Marks decahedra sequence, the (TO) + sequence and the 
truncated octahedra sequence. 

We then construct caloric curves for the lowest energy structure identified for each size 
range in the constant energy (microcanonical) ensemble using the following procedure: at 
each fixed total energy the cluster was equilibrated for at least 150000 time steps (At = 3.63 
fs) and then the kinetic energy was averaged over a further 150000 to obtain a temperature. 
An energy increment of 1.0 meV/atom was used to adjust the total energy between each 
simulation by a uniform scaling of the kinetic energy. To identify and characterise solid- 
liquid coexistence, we follow Cleveland et ali& in using the bimodality of the distribution 
of diffusion coefficients to distinguish solid and liquid atoms. This method is discussed in 
further detail in Re£i£ where it was used to characterise the coexisting solid-liquid states 
in Ag, Cu and Ni clusters. To characterise solid-solid transitions we have used common- 
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FIG. 1: Comparison of the relaxed zero temperature energies per atom of the Marks decahedra, 
truncated octahedra and the (TO) + fee closed-shell sequences for platinum clusters. Energies per 
atom are given relative to a fit (cubic in iV 1 / 3 ) to the energies of the truncated octahedra sequence. 

neighbour analysis (CNA)^ using the classification scheme given in Ref.— . 
III. RESULTS 

We have examined the thermal stability of the two closed-shell decahedra by construct- 
ing microcanonical caloric curves using the procedure described above. Fig. |2] shows the 
caloric curves for the 887-atom and 1389-atom Marks decahedra. Note that the curves are 
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qualitatively very similar with two transitions evident in each curve (indicated by the drop 
in temperature and denoted by A, B, C and D in the figure) over the range of energies 
considered. The second transition in both cases (at points B and D) is the full melting 
transition i.e. at energies above the transition point the particles are fully molten. However 
in both cases the transition prior to melting (A and C respectively) is a structural transition 
from the decahedral structure to a fee structure. From the decrease in temperature at the 
transition (indicating an increase in internal energy) we surmise that the structural transi- 
tions are not driven by energetics (i.e. the transitions are not due to a misidentification of 
the ground state structure, nor the energetic cost of phase coexistence^) . 

We will examine the transition in the 887-atom cluster in more detail (in fact we find 
similar behaviour in the 1389-atom decahedral cluster). Fig. |3] shows several snapshots of 
the 887-atom cluster, firstly at K showing the relaxed Marks decahedral structure, and 
then immediately prior to and after the transition at C in Fig |21 The bimodality of the 
distribution of mobilities indicates that prior to the transition at C the decahedral structure 
is in a coexisting solid-liquid state. The location of the liquid atoms in this coexisting 
decahedral structure (centre frame in Fig.EJ) suggests that the (lOO)-facets of the decahedron 
have melted first. The structure that appears after the transition at C (the right frame in 
Fig. |3J) also shows signs of phase coexistence, with liquid atoms covering most of the surface 
of the cluster. 

In Fig |U we show the results of a CNA structural analysis of a very long duration run 
constant energy simulation of the 887-atom cluster at E = —5.108 eV/atom (an energy just 
before point C in Fig. EJ). The cluster is initially in the decahedral structure as can be seen 
from the high ratio of bulk hep to fee atoms identified by the CNA analysis (note that in 
a relaxed 887-atom Marks decahedron at K our CNA analysis counts 360 bulk fee atoms 
and 150 bulk hep atoms). The transition from this decahedral structure to an fee structure 
appears to occurs via a melting of the decahedra between t = 2 and 3 ns where the number 
of both fee and hep atoms drops sharply. Between t = 5 ns and t = 10 ns, the number of 
fee and hep atoms remains low indicating considerable disorder. Between t = 11 and 12 ns 
there is a recrystallisation of the cluster into a fee structure as seen by the increase in fee 
atoms with no corresponding increase in hep atoms. Fig. El shows the temperature and liquid 
fraction during the simulation. The liquid fraction peaks at a value of 0.3 at times between 
the disappearance of the decahedron (t ~ 3 ns) and the appearance of the fee structure 
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FIG. 2: The caloric curves for the 887-atom (solid line) and 1389-atom Marks (dashed line) dec- 
ahedra. Both curves exhibit a transition during phase coexistence (indicated by the drops in 
temperature at points A and C respectively) prior to complete melting (which occurs at points B 
and D respectively.) 

(t ~ 12 ns). We note also that the changes in temperature correlate well with the changes 
in liquid fraction, with a minimum in the temperature occuring when the liquid fraction is 
at its highest. 
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FIG. 3: Snap-shots showing the 887-atom particle at K and then immediately prior to and after 
the transition at C in Fig. [2J The solid and liquid atoms, as identified by their mobilities, are 
shown in light and dark shades respectively. 

IV. DISCUSSION 

We have now observed several instances of solid-solid transitions that occur during solid- 
liquid phase coexistence. In one instance, a 1415-atom Ni icosahedron, a transition was 
driven by facet- dependent wetting of the melt£. In a second instance, we observed par- 
tial melting of similar sized Pd decahedra followed by recrystallisation into an icosahedral 
structure^ but this was not driven by the energetics of the of the wetting. The transition 
reported in this article is most similar to that seen in the Pd decahedra; the decrease in tem- 
perature at the transition indicates that it is not energetically driven (unlike the transition 
in Ni) and it appears to occur as a result of a melting of the decahedral structure followed 
by recrystallisation into the fee structure. In Ref.— we speculated that the transition was 
determined by recrystallisation kinetics as icosahedra are often produced by rapid quenches 
of the melt in small particles. However, the fact that the Pt particles here recrystallise into 
a fee structure rather than an icosahedron is not as easily explained by recrystallisation ki- 
netics. It seems likely that this fee structure is the stable structure at this temperature. We 
also note that the transition seen in Pd particles occured much more rapidly (recrystallisa- 
tion occured less than 0.1 ns after the decahedral structure disappeared) than the transition 
examined here (recrystallisation took 5-10 ns in the 887-atom cluster). 

It is not clear what role the phase coexistence plays in the stability of the fee structure. It 
may be that the partial melting simply facilitates the kinetics for the solid-solid transition so 
that it is observable on a timescale accessible by our simulations. On the other hand it may 
be that the fee structure is stable only in the presence of the melt. We note that solid-solid 
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FIG. 4: A plot of the number of fee (light line) and hep (bold line) atoms of the 887-atom cluster 
during a long (18 ns) constant energy simulation at E = —5.108 eV/atom. The transition from a 
decahedral structure to an fee structure appears to occurs via a melting of the decahedra (between 
t = 2 and 3 ns where the number of hep atoms drops sharply) followed by a recrystallisation into 
the fee structure (between t = 11 and 12 ns). 

structural transitions in coexisting solid-liquid clusters are not ubiquitous in clusters of this 
size (e.g. see Ref.— ) but the examples now known from molecular dynamics simulations in 
Ni, Pd and Pt clusters suggests that such transitions may be reasonably common. While 
there is no direct experimental evidence for such transitions as yet, it is possible that the 
icosahedral to decahedral transition in gold clusters observed by Koga et al^ took place in 
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FIG. 5: A plot of the instantaneous temperature (plus a time averaged trace of the temperature in 
white) and liquid fraction of the cluster during the same constant energy simulation at E = —5.108 
eV/atom as that shown in Fig. 21 The liquid fraction peaks at about 0.3 between t = 5 — 15 ns 
before dropping once more as recrystallisation occurs. 

phase coexisting particles. 

When comparing experimental determinations of nanoparticle structure with correspond- 
ing theoretical predictions, one must be aware that particles will frequently be found in 
kinetically-trapped metastable structures^. Growth through atom-by-atom aggregation in 
an inert-gas aggregation source is known to favour non-crystalline structures in certain 
circumstances^. Particles may also have undergone melting and freezing during growth, 
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as well as periods of phase coexistence 2 ^. In addition if the stable low temperature struc- 
ture becomes unstable at elevated temperatures it may be difficult to observe because of 
kinetic trapping. In the case of small Pt nanoparticles, non-crystalline structures have not 
been observed 8 -. Here we have found a process that favours crystalline structures over non- 
crystalline structures in Pt particles at elevated temperatures. It is possible that the small 
fee particles observed in Ref.— were kinetically trapped as the result of the high temperature 
instability of decahedral structures. 

V. CONCLUSION 

We have investigated the thermal stability of small decahedral platinum particles. Al- 
though these particles appear to be the stable structures at low temperatures, they undergo 
a transition to an fee structure at temperatures near the melting point. This instability may 
be why non-crystalline platinum particles have not been observed experimentally. Previous 
observations of small platinum particles^ may only have identified the metastable kinetically 
trapped fee structures. 
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